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1.0 Introduction 

In this report we examine the influence of the atmosphere 
on nadir directed signals associated with satellite altimeters . 
Frequencies at 6, 13.5, and 35 GHz are selected so as to pro- 
vide a parameter study. Uncertainties are summarized in both 
existing and proposed techniques which establish ionospheric and 
tropospheric height corrections. The error summary thus gives 
values describing the "best you can do" in height resolution 
(as dictated by atmospheric parameters) for a satellite borne 
altimeter system. The impetus to pursue this effort stems from 
a possible future satellite experiment for studying ocean surface 
currents (TOPEX) . An altimeter system aboard such a satellite 
may require mean sea height resolutions significantly better 
than 10 cm [Nagler and Brown, 1979]. 

The results presented here reflect data gleaned from the 
literature, at large, as well as from the existing body of pub- 
lished literature associated with the Seasat A Altimeter Experi- 
ment (See References - Section 6 . 0 ). 

We specifically consider here (1) the effects of pre- 
cipitation on altimeter signals (Chapter 2) (2) , range errors due 
to refractive index variations in both the clear atmosphere (con- 
vective and non-convective) and clouds (Chapter 3 ) , and ( 3 ) range 
errors introduced by the ionosphere (Chapter 4) . 

The uncertainties introduced by a two frequency altimeter 
system which continuously corrects range errors introduced by the 


Page 2 


ionosphere are examined in Chapter 4. In Chapter 2, a pre- 
liminary analysis is pursued establishing the feasibility of 
incorporating rain rate range gaces in a future satellite 
borne altimeter system. 
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2.0 Effects of Precipitation on Altimeter Signal 
2.1 Range Errors Due to Precipitation 

The range error contribution due to a rain medium 
alone may be expressed as, 

L 

AR e * f j R e ( “" 1) } dl (2>1) 

o 

where m represents the equivalent complex refractive index of 
a plane-parallel medium containing rain particles as described 
by Van De Hulst (1957), R g denotes the real part, di. the differ- 
ential path length, and L the path interval through rain (one way) . 

The complex part of m gives rise to the attenuation 
characteristics of the radio waves. In Fig. 1 are given plots 

I 

of range errors, AR^ expressed in cm per km path length through 
rain as a function of rain rate, R (expressed in mm/hr) for 6, 

13.5, and 348 GHz. The curve at 34.8 GHz was derived from the 
phase tabulations of Oguchi (1973) for the horizontal polariza- 
tion case (distorted raindrops). The curves at 6 and 13.5 GHz 
were derived from the phase results of Rogers and Olsen (1975) , 
Olsen, et al (1978) and a C.C.I.R. Document 5/23-E (1975) for 
spherical raindrops. In all case, Laws and Parsons (1943) rain- 
drop size distributions were assumed in the computation of m. 

It may be noted that similar sets of curves at other frequencies 
(1.43-300 GHz) may be derived using the equivalent refractive 
index results of Setzer (1970). In the use of such curves, a 
note of caution should be included because rain is a dispersive 
medium. Both the phase and group time delays decrease with in- 
creasing frequency, the extent of which depends on the rain. 

Rogers and Olsen (1975) point out that the theory breaks down at 
higher frequencies and the results should be used with caution. 
Examples of the dispersive nature of rain are given in Fig. 2. 






DELAY (ns/km) 
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FREQUENCY (GHz) 

Fig. 2 Time Delay Through Rain vs Frequency [Rogers, 
1975]; Tg and ip denote group and phase time 
delays, respectively. 


Ik 
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Although the results depicted in Fig. 1 are derived 
from the phase refractive indices (as opposed to the group re- 
fractive index), they are considered approximately accurate. The 
differences between the curves of Oguchi and those of Rogers and 
Olsen in Fig. 1, are due both to computational reasons as well as 
the fact that Oguchi considers distorted raindrops. 

In establishing the overall nadir rain paths , two basic 
types of rain may be considered; nonconvective (stratiform) and 
convective (thunderstorm). For stratiform rain types, the rain 
may be assumed to fall from the 0°C isotherm level. Typical 
heights of the zero degree isotherm for the band of latitudes 
ranging from 15°N to 45°N vary from 4.0 to 4.7 km during the month 
of July (Cole, et al, 1965; Cole and Kantor, 1963). Stratiform 
type rains are generally widespread (hundreds of km) and have 
rain rates of less than 20 mm/hr associated with them (C.C.I.R., 
1978). 

In convective rains , drops are often carried in the 
supercooled state to several km above the 0°C isotherm, and the 
path interval of rain may be as high as 10 km. This type rain is 
comprised of isolated intense cells or cells imbedded in rains 
covering a larger area. The horizontal extent of the intense 
regions of rain are of the order of several km (C.C.I.R., 1978) and 
many of such cells may exist within a thunderstorm region. Rain 
cell dimensions are further characterized in Section 2.4. 

We note from the above and the curves in Fig. 1, that 
for a 5 km nadir rain path and a 100 mm/hr rain (thunderstorm) , 
we obtain a 2.5 cm range error. For a stratiform type rain of 
10 mm/hr (moderate rain) over a 5 km path, we obtain a 0.5 cm 


range error. 
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Hence, we conclude that the occasional thunders tors: 
may give rise to a discernible range error whereas the more 
frequent widespread type rains may be ignored. These conclu- 
sions, however, carry with them a few caveats and these are: 

(1) the more intense rains may give rise to a considerable loss 
of signal at frequencies above 10 GHz, which may represent a 
more severe problem than that arising due to range errors (see 
Section 2.2), (2) widespread rain creates additional brightness 

temperature which interferes with the ability of auxiliary racio- 
metry measurements to amove the path error due to water vapcr 
in the troposphere (see Section 3.3). For example, the Scanning 
Multifrequency Microwave Radiometer (SMMR) on Seasat 1 was net 
always able to accurately retrieve sea surface temperatures in 
very light rains exceeding 0.5 mm/hr [Seasat Workshop Report, 

Vol . 1, p. 8-34]. 
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2.2 Attenuation Effects of Precipitation 

Raindrops both absorb and scatter radio waves giving 
rise to attenuation; absorption being the prime contribution at 
frequencies below 10 GHz. At higher frequencies, scattering be- 
comes important and at around 60 GHz and a rain rate of 5 mm/hr 
both contributions are approximately equal (Setzer, 1970). At 
the higher rain rates the scattering contribution may exceed 
that due to absorption, depending up_n the frequency. 

In Figure 3 are plotted a set of curves depicting the 
attenuation coefficients as a function of frequency and rain rates 
(C.C.I.R., 1978). These curves were obtained from the results of 
Olsen et al, (1978) asstuning a Laws and Parsons (1943) drop size 
distribution, terminal velocities of Gunn and Kinzer (1949) , and 
the index of refraction of water at 20°C of Ray (1972) . We note 
from Fig. 3 a monotonic increase of attenuation with frequency up 
to 100 GHz beyond which a leveling off and slight reduction occurs. 

In Figure 4 are plotted the attenuation coefficients 
for the specific frequencies of 6, 13.5, and 35 GHz, as a function 
of rain rate. These may be expressed by the empirical relations 
(Olsen et al, 1978), 

k = aR b (2.2) 

where k = attenuation coefficient (dB/km) at the frequency, f 
a,b = empirical frequency and drop size distribution 
dependent quantitites 


R = rain rate (mm/hr) 
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In Table 1 are given values of a and b at the frequencies given 
in Figure 4. 

Table 1 Values of a and b for k * aR^ Relation [Olsen et al, 1978] 


Frequency (GHz) 

R < 50 mm/hr 

R > 25 mm/hr 

a 

b 

a 

5 

6 

1.79xl0" 3 

1.238 

1.06xl0" 3 

1.393 

13.5 

2.48xl0' 2 

1.205 

3.46xl0' 2 

1.109 

35 

2.35X10" 1 

1.009 

3 . 37xl0 _1 

0.904 


We note from Figure 3 that at a 100 mm/hr rain rate, the 
attenuation coefficients are 0.65, 5.7, and 21.7 dB/km at 6, 13.5 
and 35 GHz, respectively. Assuming a 10 km round trip path dis- 
tance, the total respective attenuations are 6.5, 57 and 217 dB 
at 100 mm/hr rain rates. 

For Seasat A operation (f=13.5 GHz) an approximate 
14 dB margin existed before loss of track occurred (at about a 
signal to noise ratio of 6 dB) whereas optimal operation existed 
with signal to noise levels above 13 dB. For a two way 10 km 
path length, 14 dB corresponds to an attenuation coefficient of 
1.4 dB/km, which at f = 13.5 GHz, corresponds to a rain rate of 
about 30 mm/hr (Fig. 4) . 
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The passage through intense rain cells by Seasat A 
could have resulted in distorted altimeter return signals giving 
rise to subsequent misinterpretations of wave height. This is 
evidenced by the sudden apper’-ance of occasional interpreted 
wave height values up to 15 to 20 m for a distance of thirty 
kilometers or so (Lorell, 1979). An illustration of a rain cell 
encounter is depicted in Fig. 5 (Seasat Workshop Report, Vol. 1, 
p. 5-41). During the encounter we note that scattering cross 
section, o° drops about 1 dB for 4 seconds (28 km along ground 
track), the significant wave height, rises from 3 m ambient 

to a sharp peak of 6 m and drops back to ambient in approximately 
1 sec (7 km) . The sea surface height is shown to drop about 1 m 
and return to ambient in about 1 sec. 

It is apparent from the attenuation curves in Fig. 4 
and the time delay results depicted in Fig. 1, that at the rain 
rates in which the range errors become a problem for f » 13.5 GHz 
(and greater frequencies) , attenuation becomes prohibitively high, 
resulting in loss of lock or highly distorted returns . On the 
other hand, there may be marginal time delays for 6 GHz at the 
high rain rates (100 mm/hr) where the attenuation is considerably 
reduced. 


START TIME: 268:18:33:4 (IN SECS 23135585.0742) 



Response of Seasat Altimeter to a Rain Cell (Seasat 
Gulf of Alaska Workshop Report, Vol. 1, April 1979) 
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2.3 Measurement of Earth's Rain Intensity with a Satellite 

Altimeter 

We examine here the power received at the input of a 
Seasat type altimeter due to rain scattering at frequencies of 
6, 13.5, and 35 GHz. The measurement of rain and the corre- 
sponding intensity (mm/hr) could serve multifold purposes; namely, 

(1) Determination of the existence of rain rate 
levels could explain deterioration of altimeter 
performance (e.g., increased noise, and higher 
uncertainties in mean sea level, significant 
wave height, and backscatter coefficient of 
the ocean surface. 

(2) Explain poor radiometer performance and estab- 
lish regions where the atmospheric corrections 
using radiometry should not be used. 

(3) Provide rain rate intensities and vertical pro- 
files on a world-wide scale. 

2.3.1 Detectable Power from Rain Scatter for Seasat Type 
Altimeters 

The power scattered from a pulse volume of raindrops 
is given by (Goldhirsh, 1979) 



o 


(2.3) 


where 
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c 



T 


X 

G 


A0 e ,A0 h 


L 


r 


W k = 


velocity of light (3 x 10 m/s), 
transmitted power (watts) , 
pulsewidth (sec) , 
wavelength (m) , 
antenna gain, 

vertical and horizontal beamwidths (rad) , 
transmitted loss factor (<,1.0) (losses from 
transmitter power measurement point to gain 
measurement point) , 

receiver loss factor (11.0) (losses from 
receiver calibration point to antenna gain 
measurement point) , 
rain reflectivity (m~^) , 
range (m) , 

attenuation coefficients due to atmospheric 
gas, precipitation, and clouds, respectively 
(dB/km) , 


The factor f (B) represents the additional loss factor due to 
the radar receiver frequency response. For a matched filter 
linear receiver (bandwidth equal approximately to reciprocal 
of pulsewidth), f(B) has been estimated previously as equal to 
0.66 (-1.8 dB) by Nathanson and Smith (1972). More recently, 
f(B) has been calculated by Doviak and Zrnic (1979), for a 
Gaussian shape receiver response and rectangular transmitted 
pulse, and is given by 
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f(B) 


coth 


ttBt 

2 V ln 2 



-1 

ttBt \ 

V In 2 ) 


(2.4) 


This results in a -2.3 dB loss for the "matched filter" case, 
where B is the receiver bandwidth (Hz) . 


At the frequencies considered and for nadir directed 

rf energy, we may simplify (2.3) by allowing k as 0 and k ^ 0 

6 ^ 

(Goldhirsh 1979) . We also assume the receiver bandwidth is 
sufficiently large such that f(B)*l. For the case in which 
Rayleigh theory is applicable (e.g., f < 6 GHz; the raindrop cir- 
cumference to wavelength ratio <<c l) , Rayleigh scattering theory is 
applicable and 


n - (tt 5 / A 4 ) K q 2 Z (nf 1 ) 

where 


(2.5) 


D 

max 

Z - / N(D)D 6 dD (m 3 ) (2.6) 

D min 

The parameter Z is called the reflectivity factor and 

is defined by (2.6), where N(D)dD is the number of drops whose 

diameters are between D and D + dD per unit volume (drop size 

distribution DSD), and D max » D m i n are the maximum and minimum 

diameters of the drops sampled in the pulse volume, respectively. 

2 

The parameters |K q | is dependent on the complex refractive index, 
and has been tabulated for various frequencies and temperatures 
by Gunn and East (1954) and Ray (1972) . In the frequency range 
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6 Co 35 GHz and for temperatures ranging from 0° to 20°C, 

|K 0.9. At 13.5 GHz Rayleigh theory may not be applicable 
and a 407. underestimation in reflectivity may result if (2.5) 
is used (Goldhirsh, 1977). However, within the uncertainties 
considered here we shall assume (2.5) is valid at both 6 and 
13.5 GHz. At f-35 GHz (Goldhirsh, 1977; p. 93) 

n 4.43 x IQ" 7 Z 0 * 745 ; f-35 GHz (nf 1 ) (2.7) 

6 3 

where Z assumes the more convenient units in (2.7) of (mm) /m . 

Assuming the altimeter parameters at the three fre- 
quencies considered conform to those on Seasat 1; namely, 

P - 2 kw 

t - 3.2 sec 

r - 800 km (2.8) 

A0 «A0 l - 1.6° 

e h 

G - 40.6 dB 

and considering the aforementioned simplications , (2.3) becomes, 
P r (dBm) - A + B log 1Q R - 2 aR b L (2.9) 

where the Z's in (2.5) and (2.7) have been related to R by , 

1.6 


Z - 200 R 


( 2 . 10 ) 
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a relationship which has been found to approximately conform 
to the Laws and Parsons (1943) drop size distribution (approxi- 
mate stratiform type rain). Also kp in (3.1) has been related 
to R via the relation (2.2). That is, 

k p - aR b (2.11) 

and L in (2.9) is the nadir directed path length (one way) through 
rain of rain rate, R, assumed uniform in the vertical extent. 

In Table 2 are given values of a and b (Olsen et al, 1978) 
and A and B for 6, 13.5 and 35 GHz used in (2.9). 

Table 2 Values of A, B, a and b in the Relationship (2.9) 


Frequency (GHz) 

A 

B 

a 

b 

6 

-120 

16 

1 . 79xl0 -3 

1.238 

13.5 

-113 

16 

2 . 48xl0 -2 

1.205 

35 

-101.2 

11.9 

2 . 35xl0 -1 

1.009 


In Figs. 6 and 7 are plotted the altimeter received 
power as a function of rain rate as given by (2.9) where > rig. 6 
we compare the 6 and 13.5 GHz cases and in Fig. 7 the 12 •' ..*d 35 
GHz cases. The straight lines correspond to the case of negligible 
attenuation and are associated with measurements at the rain tops 
(i.e., L -0) and the curved lines depict the power received from 
a pulse volume at the earth’s surface where L-5 km. 
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The bending of the power curves in Figs. 6 and 7 are 
due to two opposing factors in the power relationship as the 
rain rate increases. At larger rain rates the backscatter re- 
flectivity increases tending to increase the power. However, the 
path attenuation also increases tending to reduce the power. At 
rain rates larger than those defined at the curve maxima, the at- 
tenuation term dominates over the backscatter reflectivity term 
and a net reduction in power results . 

A nominal value of noise for the satellite altimeter 
is given by, 


N = k T B F > - 109 dBm (2.12) 

where 

k = 1.38 x lO' 23 Joules/°K (Boltzmann's Constant) 

T = 310°K 

F = 10 (Noise Figure of Seasat 1) 

B ^ i = 3.125 x 10 5 Hz 

We thus note from Figures 6 and 7 that the power scattered from 
the rain regions may be well above the receiver noise level. For 
example, a dual wavelength system at 6 and 13.5 GHz would enable 
the measurement of rain rates above 2 mm/hr (Fig. 6). On the 
other hand a 35 GHz system would enable rain measurements between 
0.2 and 8 mm/hr (Fig. 7). 

Improvement of the noise figure in future systems (e.g. , 
F -*■ 5 dB) should result in further enhancement of the signal and 
enable the determination of a larger range of rain rates . 


2.4 Rain Cell Dimensions 
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Considerable ambiguity exists as to the diameter of 
rain regions. Many investigators have characterized cell di- 
mensions (C.C.I.R. , 1978). however these values are very much 
dependent on the definition. In actuality, both stratiform type 
and thunderstorm rains often contain many embedded regions of 
intense rainfall. Hence, peaks and valleys of rainfall intensity 
often exist along the horizontal direction, which makes it dif- 
ficult to relate a rainrate to a rainfall dimension. Notwith- 
standing this type of ambiguity, we depict two rainfall models 
in Fig. 8 (taken from Skolnik, 1974) giving the storm diameter 
as a function of rain rate; namely that of Edgar et al (1973) and 
Nathanson (1969) . We note that 20 km rainfall diameters and larger 
coexist with rainfalls of 10 to 15 mm/hr and less. (This compares 
with a footprint of 22 km for a 1.6° beamwidth at 800 km range.) 
Rainfalls having larger rainrates are confined within smaller dia- 
meters and may not fill the pulse volume. A reduced value of 
scattered power may result for such cases due to averaging within 
the pulse volume. 

Median values of rainfall rates as a function of 
height have been described by Goldhirsh and Katz (1979) based on 
a large sampling of radar measurements (Konrad; 1978) . These 
median values, which are plotted as a function of ground category 
level of reflectivity factor in Fig. 9, show that within +2.5 
dBz (d3z = 10 Lc$ Z [mm /m j) , the rainfall rat is are relatively 
uniform up to approximately the zero degree isotherm, above which 
they fall off rapidly. The fall-off of reflectivity occurs at 
the higher altitudes with increased ground values of Z probably 
due to convective updrafts carrying the precipitation to higher 
altitudes. 


STORM DIAMETER 


DIAMETER OF RAINSTORM Page 

VERSUS RAINFALL RATE 



RAINFALL RATE - mm/hr 


Fig. 8 Diameter of Rainstorm Vs Rainfall Rate 



ALTITUDE (km) 
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RAIN RATE, mm/hr 



Fig. 9 Median Reflectivity Factor Profiles for Given 

Ground Categories as Measured at Wallops Island, 
Virginia [Goldhirsh and Katz; 1979] 
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3.0 Height Corrections Due to Refractive Index Variations in 
the Troposphere 

3.1 Clear Air Regions - Formulation 



N - 10 6 (n-1) 


(3.2) 


where N may be approximated by [Smith and Weintraub, 1953] 

N » 77.6 | + 3.73 x 10 5 (3.3) 

where , 

P = total atmospheric pressure [mbar] 
e - partial pressure of water vapor [mbar] 

T - atmospheric temperature [°K] 

The constants in (3.3) are considered valid to 0.5% 
in N for frequencies up to 30 GHz and normally encountered ranges 
of pressure, temperature, and humidity and with the exception of 
the oxygen constituent up to 72 GHz [Bean and Dutton, 1968]. Sub- 
stitution of (3.3) into (3.1) results in, 
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H 


H 


AH - 77 


6 f\ 


dh + 3.73xl0 5 f K dh 

T J 

o o 


or 


AH - ( AH )pRY + ^ AH ^WET 


(3.4) 


(3.5) 


where the dry and wet parts refer to the integral expressions 
of the first and second terms in (3.4). 

3.2 Dry and Wet Terms as * Function of Surface Parameters 

The height correction as given by (J.4) requires a 
knowledge of P, T, and e along the propagation path. Although 
those quantities may be obtained using radiosonde data, they are, 
in general, unknown, and simplifications of the dry and wet parts 
of (3.4) in terms of the surface parameters are desirable. The 
Seasat experimenters [Seasat Users Handbook, 1979] employed the 
Saastamoinem model [1972], given by, 


(AH) WET " 2 ‘ 277 x 10 ‘ 3 C-05 + ±f^] e g [m] (3.6) 

(AH) DRY " C2 . 277 x 10" 3 - 1.11 x 10~ 5 cos*]P g [m] (3.7) 


where the subscript s denotes the surface parameters and <P is the 

latitude (T is in units of °K and e„ and P o are in mbar) . 
s s s 
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In the determination of the wet term, (3.6), the following 
assumptions are made: 

(1) The temperature reduces with altitude at a constant lapse 
rate given by, 

T - T + 6 (h - h ) (3.8) 

8 8 

HT 

where ^ ■ constant temperature lapse rate 

h - surface altitude above sea level. 

8 

(2) The partial pressure of water vapor varies with altitude 
according to the relation, 



where k ■ ^ (3.10) 

where g - gravity constant 

R * gas constant of dry air 

In the absence of SMMR data, Seasat experimenters used 
interpolated surface values obtained primarily from the Fleet 
Numerical Weather Center (FNWC) . 

3.3 Seasat Scanning Multi-Frequency Microwave Radiometer 
[SMMR] Measurements of the Wet Term 

The Seasat radiometer system operated at frequencies 
of 6.6, 10.7, 18, 21 and 37 GHz at vertical and horizontal 
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polarizations. Using these measurements the integrated columnar 
water vapor mass density, W Cgm/cnr] may be measured. That is, 

H 

W - J P y dh Cgm/cm^l (3.4) 

o 

3 

where P v is the water vapor density in gm/cm . 

Assuming the ideal gas law, an atmosphere of 
constant density, and a temperature profile with a constant 
lapse rate, the wet height correction term has been shown to 
be given by [Lipes, 1979; Seasat User's Handbook, 1979], 


Ah wet 


* W ^(t)" k 3 (I)(t^) ln ( X 



(3.12) 


where 


W 

R 

M 

T 


k 2 ,k 3 


the SMMR measured integrated water vapor density 
ideal gas constant 
molecular weight of water 
surface temperature (air) 

temperature difference between the surface and the 
scale height 

constants related to the wet terms of the 


refractive index 
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3.4 Comparisons of Seasat Range Corrections with those 

from Radiosonde Data 

In order to arrive at measures of uncertainties for 
the altimeter height corrections, we compare here actual 
measurements taken during Seasat operation with those obtained 
using radiosonde data. For one set of data we examine four 
passes over Bermuda (passes 1117, 1160, 1246, 1375) and for the 
second set, we examine the results of four additional passes 
(1163, 1212, 1292, 1298) taken aboard the Oceanographer Vessel 
near the Gulf of Alaska taken as part of the GOASEX (Gulf Of 
Alaska Surface Experiment) . The original data examined here were 
extracted from the Seasat Gulf of Alaska Workshop Report, Vol. 1, 
April 1979. Radiosondes were taken during the eight passes to 
provide "truth measurements". 

Two algorithms were used for converting the SMMR data 
to altimeter path length corrections and these are known as the 
Wilheit and the Wer.tz algorithms. In Fig. 10, we compare the 
results of Wentz and Wilheit with those obtained using radiosonde 
data for the 8 passes mentioned above [Seasat Gulf of Alaska 
Workshop Report; 1979] . Those radiosonde points labeled "est" 
correspond to cases in which the measured profiles terminated 
prematurely and estimated corrections were employed. 

In Tables 3 and 4 we summarize the results for the 
Bermuda passes (Table 3) and the combined Bermuda-Gulf of Alaska 
passes (Table 4) . 



Fig. 10 Comparison Between Altimeter Path Length Correction 
Algorithms and Radiosonde Data (Seasat Gulf of 
Alaska Workshop Report; 1979) 





COMPARISON OF RADIOSONDE WITH MODEL AND 
SMMR HEIGHT CORRECTIONS 
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The first two columns of Table 3 represent the rms 
and % rms differences and the last two columns the maximum 
and 7. maximum differences. In the first three rows the dif- 
ferences are taken relative to the radiosonde truth corrections, 
and in the last row the theoretical models for the dry term are 
compared using the FNWC and surface measurements; the latter 
representing the assumed "truth" values. In comparing the 
SMMR measurements, the average values deduced from the Wentz and 
Wilheit algorithms were used. In Table 4 the uncertainties in the 
Wentz and Wilheit algorithms are individually compared with the 
8 radiosonde passes. 

In Fig. 11 the height corrections using the Wentz and 
Wilheit algorithms are compared with each other and with the 
theoretical model using FNWC estimates for Revolution 1375. It 
is interesting to note that a 2 cm rms difference exists when 
the Wilheit correction is compared to the reference Wentz case 
(as an example) and an approximate 6 cm rms difference exists 
upon comparing the theoretical model using FNWC surface estimates 
with the average of the Wentz an 3 lheit cases. 

Based upon the results in Tables 3 and 4 and Fig. 11, 
we conclude that the radiometer measurements give rise to a nominal 
rms height correction uncertainty of 2.5 cm and the theoretical 
models uncertainties for the wet and dry terms using FNWC inter- 
polated surface values are approximately 6 cm and 2 cm, respectively. 



WO NOIJ.naiHJ.NOO 13M 


DAY 274, REV 1375, LAT. LONG. DEG 

Fig. 11 Comparison of SMMR Algorithms and Tropospheric 
Correction Model (FNWC Data) - (Seasat Gulf of 
Alaska Workshop Report ; 19 79) . 
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3.5 Height Correction Fluctuations in the Clear Atmosphere 
and Due to Clouds 

Clouds 

In Figs. 12 and 13 we show two cases of refractive index 
variations through clouds [Campen et al. 1961] obtained with a 
refractometer located aboard an aircraft. The indicated refrac- 
tive index values represent the averages of approximately 30 
passes through the clouds at the various indicated heights for 
Boston, Mass. (Fig. 12) and Tucson, Arizona (Fig. 13). In the 
former case, the cloud extends from approximately 1.8 to 3.5 km 
and has a maximum horizontal dimension exceeding 1 km. Relative 
to regions outside the cloud, refractive index changes as high 
as 25 N units appear over a 1 km horizontal extent. For the case 
of Fig. 12, integrating N as a function of altitude results in a 
localized height of 2.3 cm. 

The cloud in Fig. 13 extends from 3.8 to 5.8 km with a 
maximum horizontal extent of approximately 2.5 km and a maximum 
N (relative to the outside of cloud) of 10. integration of N over 
the altitude extent of the cloud for this case results in a 1.8 
cm fluctuation of the height. 

We thus note that cumulus type clouds may, in fact, 
result in fluctuations of AH of a few cm. It is interesting to 
note that bands of such clouds may develop in the form of cloud 
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streets as depicted in Fig. 14 which is a photo as seen from 
Appolo 6 C Kuettner , 1971] . The length of the streets may extend 
over ICO km with a spacing of 2 to 2.5 km. 

Clear Air Convective Regions 

Similar type variations in refractive index may occur 
in the clear air ac the boundaries of moist and dry air. An 
example of aircraft measurements of the refractive index varia- 
tions at Wallops Island. Virginia is given in Fig. 15 
I Konrad. 1970a] . Intensity modulated displays (RHI's - vertical 
scans, and PPI's - azimuthal scans) of clear air convective cells 
taken at Wallops Island, Virginia, are given ir. Fig. 16 and 17 
I Konrad, 1970b] . From figures of these type, it may be demon- 
strated that £H variations and spatial structures are similar 
to those for cumulus clouds. A typical model for Wallops Island, 
Virginia of the £5i variation for a convective field is given in 
Fig. 18. 

In Tacle 5 is given a summary of the scale dimensions 
for cloud and clear air convective regions. If the footprint 
of an altimeter is within the scale dimensions indicated here, 
rapid fluctuations of the height correction of a few cm may resul 
other vise smoothing occurs. 


I 
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Fig. 14 
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Fig. 16 Radar RHI in Clear Air Convection at Wallops Island, 
Virginia [Konrad, 1970] 
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Fig. 17 Radar PPI in Clear Air Convection at Wallops Island, 
Virginia [Konrad, 1970] 
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TABLE 5 

SUMMARY OF NOMINAL SCALE DIMENSIONS FOR THE TROPOSPHERE 


CUMULUS .CLOUD 

1) MEAN DIAMETERS 

0,1 TO 6 

2) CLOUD STREETS 

LENGTH 

SPACING 

SPACING-TO-HEIGHT RATIO 


(LOPEZ; 1977) 

i 

(KUETTNER; 1971) 

20 TO 500 km 

2 TO 8 km 
2 TO 4 


clear a i r mm 

1) CELL DIAMETERS 

2) CELL FIELDS 


APPROXIMATELY SAME AS CUMULUS 
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4.0 Range Errors Introduced by the Ionosphere 

4.1 Range Error Correction 

The excess group delay over a nadir propagation path, 
H, through the Ionosphere relative to that same free space path 
is, 


n n 

4T - ^ y*( n-l) ds - ^y*4n ds 


(4.1) 


where c - velocity of light 

n ■ group refractive index 

For propagation at frequencies considered here greater 
than 1 GHz), An is given by [Lawrence et al . , 1964], 


An ^ 40.3 (4.2) 


3 

where, N - electron density (electrons/m ) 
f ■ frequency in Hertz 


Combining (4.2) and (4.1), the excess time delay in- 
troduced by the ionosphere over a range, R, is 


AT - 1.34 x 10 -7 [sec] (4.3) 

f z 


where N^, is the ionospheric electron content defined by, 


The excess range correction introduced by the ionosphere, Ah^, 
is, 

N t 

AH. - c AT - 40.3 -i- Cm] (4.5) 

1 f 2 

In Fig. 19, we plot the correction, Ah^ as a function 

of values of N ( ranging from those nominal (daytime) values at the 

16 -2 

minimum of the solar cycle (N t 20 x 10 m ) to near the maximum 
(N t 100 x 10^ m”^) for the frequencies of 6, 13.5 and 35 GHz. We 
note at the indicated minimum solar cycle, the respective range 
errors are 22.2, 4.4, and 0.7 cm and at the maximum solar cycle they 
are 112, 22.1, and 3.3 cm. These values establish nominal range 
corrections to be accounted for in any altimeter system. 

4.2 Ionospheric Scale Size 

Utilizing combinations of geostationary and orbiting 
satellite data, real time contours of electron content maps have 
been constructed by Davies [1978] and Davies et al [1977] (for 
near solar minimum), an example of which is depicted in Fig. 20. 

The indicated numbers when multiplied by 10^ represent the 
electron content values. We note that at latitudes between 20° 
and 30° and during mid- and late afternoon hours, relatively 
high gradients in the electron content may result. These trans- 
late to values as high as 2 cm/ 100 km at 13.5 GHz for the example 
cited. 
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In addition to scale changes due to diurnal variations, 

ionospheric disturbances may result in as much as a 20% variation 

of the TEC over a 100 km distance in the mid-latitudes (+ 20°)*. 

During a solar maximum this gradient could translate into 4 cm/ 100 

16 2 

km for an electron content of 100 x 10 electrons/m at 13.5 GHz. 


4.3 Faraday Rotation Technique 

A technique employed for the Seasat A altimeter experi- 
ment to account for ionospheric group delay was to measure N T 
using the Faraday Rotation Technique. The amount of Faraday 
Rotation, over a nadir range, H, is given by [Bums and Fremouw, 
1970], 

H 


= 2.36 x 10‘ 


/ 


B cos6 N ds [radians] (4.6) 


where 0 is the angle between the wave normal and the earth's mag- 
netic field (expressed above in gammas) , N is the electron density 

3 

in electrons/m , and H is in meters. 

A disadvantage of the Faraday Rotation technique is 
the fact that the quantity 

H 

B cos© N ds (4.7) 


/ 


is measured rather than 


^Private Communication cn 4/18/80 with H. Soicher, U.S. Army 
Communi cat ions R and D Command, Ft. Monmouth, New Jersey. The 
above comment is based on Dr. Soicher 's measurements of electron 
content . 


Page 50 


H 

J N ds (4.8) 

0 

In arriving at (4.8), the term B cos0 in (4.7) is usually re- 
moved from the integral sign and replaced by it 3 value at a 
"mean" ionospheric altitude (e.g.,—400 km (Klobuchar ; 1978) ) . 

Since B decreases inversely with the cube of the geocentric dis- 
tance, and the electron density decreases exponentially with alti- 
tude above F 2 max (i^300 km), the integral is heavily weighted near 
the earth and is considered to provide electron content values for 
altitudes below ^1200 km [Soicher 1975]. 

For the Seasat altimeter operation, Faraday rotation 
measurements were made at Goldstone, California and Armidale, 
Australia using the ATS-1 and ATS-2 geosynchronous satellite 
beacons. Total electron content measurements, N T , were deduced 
at these stations and extrapolated to sub-satellite points on 
earth from an algorithm [Seasat User's Handbook, 1979; VJu, 1977] 
having the functional dependence given by, 


Nip f {Np(min) , Nip(max) , x» b } 


(4.9) 


where Np(min), Np(max) = minimum and maximum electron contents 
as measured at the stations 


= g eoma 8 ne bic latitude 
X = sun angle 
t = local time 
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The accuracy of N T prediction by extrapolating 
Faraday measurements at a single location to arbitrary sub-satellite 
points may be as high as 507o [Klohuchar; 1978], the actual level of 
uncertainty depends upon the year within the solar cycle, location, 
season, as well as the time of day. It may be noted that the most 
recent solar minimum occurred around 1974 [Oh; 1974]. The standard 
deviation of values of electron content at any one location based 
on daily measurements over any one month was found to be nominally, 
25% [Klobuchar , 1978] . 

In an accuracy comparison between the Seasat algorithm 
and the model used by the Naval Service Weapon Center (NSWC) , differ- 
ences as much as 5 cm were noted for a Seasat pass on 9/13/78 
[Tapley et al; 1979]. The latter model used daily measurements 
of the sunspot number (as measured at Boulder, Colorado) to obtain 
the electron content. 

4.4 Two Altimeter Method Correction 

A method by which continuous sub-orbital point 
ionospheric corrections may be obtained is described here. This 
method employs two altimeter systems at frequencies fy and fy 
(e.g., fy = 6 GHz and fy = 13.5 GHz; L and U denoting lower and 
upper frequencies). In Fig. 21 we depict the corresponding 
idealized transmitted pulses and the respective mean return wave 
forms at frequencies fy and fy. The times Aty and Aty represent 
the respective delays between the leading edge points of the trans- 
mitted pulses and the indicated leading edges of the mean return 
waveforms [Townsend, 1980]. 


1 
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Fig. 21 Two Frequency Altimeter Method for Ionospheric 
Height Correction. Subscripts L and U denote 
Lower and Upper Frequencies, respectively. 


l 
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These delays are related to the true range, H^,, In 
the following manner. 


- 4H 1U 

(4.10) 

- iH iI. 

(4.11) 


where Hy and Hy are the apparent measured ranges at fy and fy 
given by, 


H 


U 


«L 


(C Aty) 


(C Aty) 


(4.12) 


and AHyy and AH^y are the ionospheric height corrections given 
by (4.5) at fy and fy, respectively. We note that, 


( iH lL| / f U^ 

\^l \*l) 


K 


(4.13) 


Injecting (A.. 13) into (4.10) and (4.11), and eliminating AH^y and 
AHiy, we obtain 


h t = (k7t) H U " (iTl) ^ 


(4.14) 
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Allowing and cr^ to represent the system standard 
deviations associated with the Hy and measurements , respec- 
tively, it may be demonstrated that in subtracting out the 
Ionospheric contributions and arriving at (4.14), the resultant 
uncertainty in H^, becomes, 


a 


i 


M 


'U 


(4.15) 


where , 


M - 


ITT 



(4.16) 


In (4.15) and (4.16), M derotes a multiplication factor introduced 
by the additional lower frequency altimeter system. We note that 
as the frequency separation increases, K assumes large values 
and M — > 1. In Fig. 22 are plotted values of M as a function 
f L where it is assumed that fy - 13.5 GHz. The different curves 
correspond to various assumed uncertainty levels in tha measure- 
ment of H^. We note that at ■ 6 GHz and ^ = °y = 2 cm 
(Seasat case for °y) , 


* 2.5 cm (4.17) 

The two frequency altimeter system thus eliminates the range 
uncertainty due to the ionosphere although an increase of error 
of 0.5 cm (e.g., from 2 to 2.5 cm) does occur as compared 
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to that of the single frequency case (for the example considered) . 

An important byproduct of the two altimeter system Is 
that it would provide highly accurate measurements of continuous 
swaths of the ionospheric electron content. Such measurements 
would be capable of "filling in" large gaps in regions over the 
oceans where few measurements exist and hence, establish measured 
worldwide scale size information. It would thus provide to the 
engineer information on range error variability necessary for 
the design of future high precision altimeter missions. 
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5.0 Summary and Conclusion 

In Table 6, we summarize the nominal uncorrected 
atmospheric height errors for the troposphere and Ionosphere 
In the absence of rain effects. We note the dry tropospheric 
term due to all the atmospheric constituents with the exception 
of water vapor, Is significantly larger than the wet term, due 
only to water vapor. However, the variability of the dry term is 
significantly less than that of the wet term and the uncertainty 
in predictability is significantly smaller. The ionospheric cor- 
rection depends on the time of day, location, season, and year 
with the solar cycle. Nominal day time mid- latitude values are 
shown at the three frequencies where "MIN" and "MAX" denote electron 

contents of N T ■ 20 x 10*^ electrons/m^ (solar min) and 100 x 10^ 

2 

electrons/m (solar max) . 

In Table 7, we summarize the height uncertainties (rms) 
after the implementation of atmospheric corrections assuming the 
indicated modes of operation as described by column (1) . For 
example, Row A denotes a satellite containing a two frequency 
altimeter (13.5 and 6.0 GHz), and a radiometer system. Column 
(2) indicates "rain" or "no rain", the former case generally 
disallowing radiometer measurements of the wet term. Column (3) 
corresponds to uncertainties in removing the range errors for 
the wet term in the troposphere when a radiometer is used, and 
columns (4) and (5) apply when theoretical models are incorporated. 
Column (6) gives the additional uncertainty in the removal of the 
ionospheric correction when a two frequency altimeter is employed 


ARY OF NOMINAL ATMOSPHERIC UNCORRECTED HEIGHT 
ERRORS FOR THE TROPOSPHERE AND IONOSPHERE 




















Table 7 

nominal height uncertainties after atmospheric corrections are implemented 
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and (7) when employing the Seasat ionospheric algorithm. Column 
(8) represents the total uncertainty (square root of the sum of 
squares of individual errors) . 

In comparing the various total uncertainties in Table 7 
(rows A and C) , we note that the employment of a two frequency alti- 
meter system results in an nominal net atmospheric precision of 3 cm 
as compared to 6 cm assuming the Seasat system. Since the ionospheric 
correction is continuously and systematically implemented along track 
with the two frequency method, significantly larger peak uncer- 
tainties (e.g. , 5 cm) are expected to be removed. A byproduct 
result of such a system would be to provide (for the first time) 
a world wide continuous swath of the electron content (integrated 
columnar electron density) . Aside from providing pertinent data to 
ionospheric physicists, such data would be extremely useful in 
establishing scale variabilities of ionospheric, electron content; 
a useful input in the design of future altimeter missions. 

The feasibility of gating rain reflectivity with the 
altimeter system is established in Chapter 2. Such a capability 
would provide the experimenter with an accurate estimate of rain 
rates (mm/hr) at various altitudes. Such results could be employed 
in culling of the data, providing ancillary information to the al- 
timeter measurements. A continuous track of rain rates values could 
also provide the meteorologist with quantitative information on a 
world wide scale over the oceans where relative little rain measure- 


ments exist. 
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In Incorporating a radiometer system for removing 
the wet contribution, the radiometer antenna beams should gen- 
erally be coincident with that of the altimeter Bystem. This 
would tend to Insure that the height corrections apply to the 
same column of the troposphere unlike the SMMR on Seasat. (The 
SMMR had variable footprints of 21 km at 37 GHz to 121 km at 6.6 
GHz and the main beam was offset from nadir by approximately 50° 
[Dunne, 1978].) 

Also suggested in a future high precision altimeter 
mission is a nadir pointing coincident infrared radiometer for the 
detection of clouds (e.g., at 10 microns). This device would prove 
highly useful in establishing correlation between possible anomalous 
height variations and the presence of cloud [Grody, 1976]. 


Page 62 


6.0 References 

Bean, B.R. and E. J. Dutton, [1978], "Radio Meteorology", Dover 
Publications, New York. 

Burns, A. A. and E.J. Fremouw, [1970], "A Real Time Correction 

Technique for Trans- ionospheric Ranging Error", IEEE Trans, 
on Anten. and Prop., Vol. AP-18, No. 6, November, pp. 785-790. 

Campen, C.F. ,Jr. , R.M. Cunningham, V.G. Plank [1961], "Electro- 
magnetic Wave Propagation in the Lower Atmosphere" , Chapter 
13 of Handbook of Geophysics, MacMillan Co., New York, pp. 

13-5 to 13-11. 

C.C.I.R International Radio Consultative Committee, [1978], 

Recommendations and Reports of the C.C.I.R., XIVth Plenary 
Assembly, Kyoto, Vol. V, Propagation in Non-Ionized Media. 

C.C.I.R. Document 5/23-E, Canada, 22 October [1975], Study Program 
5C-2/5, "Propagation Data Required for Space Telecommunication 
Systems" . 

Cole, A.E., A. Court and A.J. Kantor, [1965], "Model Atmospheres" 
in Handbook of Geophysics and Space Environments, edited by 
S.L. Valley, pp. 2-11 - 2-13, McGraw-Hill, New York. 

Cole, A.E. and A.J. Kantor [1963], Air Force Interim Supplemental 
"Atmospheres to 90 Kilometers", AF Surv. Geophys., No. 153, 

Air Force Cambridge Research Laboratories, Hanscom AFB, Mass. 

Cunningham, R.M. [1964], "Meteorological Aspects of Range and 
Range Rate Error", AFCRL Report ESD-TO-4-64-103 . 

Davies, K. [1978], "Forecasting ami Prediction of Ionospheric 

Parameters", AGARD Lecture Series No. 93, Recent Advances in 
Radio and Optical Propagation for Modern Communications, 
Navigation and Detection Systems, pp . 6-1 to 6-29. 

Davies, K. , G.K. Hartmann and R. Leitinger [1977], "A Comparison 
of Several Methods of Estimating the Columnar Electron Con- 
tent of the Plasmasphere" , J. Atmos. Terr. Res., Vol. 39, 
pp. 571-580. 


Page 63 


References (Contd.) 

Doviak, R.J. and D. Zrnic [1979], "Receiver Bandwidth Effect on 
Reflectivity and Doppler Velocity Estimates", J. Appl. 
Meteorol., Vol. 18, No 1. pp. 70-76, January. 

Dunne, J.A. , [1978], "The Experimental Oceanographic Satellite 

Seasat-A", Boundary- Layer Meteorology, Vol. 13, pp. 393-404. 

Edgar, A.K. , E.J. Dodsworth and M.P. Warden [1973], "The Design 
of a Modern Surveillance Radar", International Conference on 
Radar - Present and Future, 23-25 October, IEE Conference 
Publication No. 105, pp. 8-13. 

Essen, L. and K.D. Froome [1951], "The Refractive Indices and 
Dielectric Constants of. Air and Its Principal Constituents 
at 24000 Mc/s", Prec. Phys . Soc. 'London, England), Vol. 64, 
October, pp. 862-875. 

Goldhirsh, J., [1979], "A Review on the Application of the Non- 

Attenuating Frequency Radars for Estimating Rain Attenuation 
and Space Diversity Performance", IEEF Trans, on Geoscience 
Electronics, Vol. GE-17, No. 4, October, pp. 218-239. 

Goldhirsh, J. and I. Katz, [1979], "Useful Experimental Results 
for Earth-Satellite Rain Attenuation Modeling", IEEE Trans, 
on At;ten. and Prop., Vol. AP-27, pp. 413-415. 

Goldhirsh, J., [1977], "Analysis of Single and Multiple Wavelength 
Radar Systems for Estimating Slant Path Attenuation Through 
Rain and Cloud at Frequencies Greater than 10 GHz", APL/JHU 
Technical Report S1R77U-036, October. 

Grody, N.C., [1976], "Remote Sensing of Atmospheric Water Content 
from Satellites Using Microwave Radiometry", Trans, on Anten. 
and Prop., Vol. AP-24, No. 2, March, pp . 155-162. 

Gunn, K.L.S. and T.W.R. East [1954], "The Microwave Properties of 
Precipitation Particles", Quart. J. Roy. Meteorol. Soc., Vol. 
80, pp. 522-545. 


rage 04 


References (Contd.) 

Gunn, K.L.S. and G.D. Kinzer, [1949], "The Terminal Velocity of 
Fall of Water Droplets in Stagnant Air", J. Meteorology, 

Vol. 6, pp. 243-248. 

Klobuchar, J.A., [1978], "Ionospheric Effects on Satellite Navi- 
gation and Air Traffic Control Systems", AGARD Lecture 
Series No. 93, Recent Advances in Radio and Optical Propag. 
for Modern Communications, Navigation and Detection Systems, 
pp. 7-1 to 7-17. 

Konrad, T.G. , [1978], "Statistical Mou^ls of Summer Rainshcwers 
Derived from Fine-Scale Radar Observations", J. Appl. 
Meteorol. , Vol. 17, No. 2, pp. 172-188, February. 

Konrad, T.G. and F.L. Robison, [1973], "Development and Charac- 
teristics of Free Convection in the Clear Air as Seen by 
Radar and Aircraft", J. of Appl. Meteorol., Vol. 12, No. 8, 
December, pp. 1284-1294. 

Konrad, T.G. , [1970a], "The Dynamics of the Convective Process 
in the Clear Air as Seen by Radar", 14th Rader Meteorology 
Conference, November 17-20, Tucson, Arizona, American 
Meteorological Society, pp. 57-67. 

Konrad, T.G. , [1970b], "The Dynamics of the Convective Process 

in Clear Air as Seen by Radar", J. oi atmospheric Sciences, 
Vol. 27, No. 8, November, pp. 1138-1147. 

Kuettner, J.P., [1971], "Cloud Bands in the Earth's Atmosphere - 
Observation and Theory", Tellus, Vol. 23, pp. 404-425. 

Lawrence, R.£., C.G. Little and H.J.A. Chivers , [1964], "A 
Survery cn Ionospheric Effects Upon Earth-Space Radio 
Propagation", Proc. IEEE, Vol. 52, pp. 4-27, January. 

Laws, J.O. and D.A. Parsons, [1943], "The Relation of Raindrop 
Size to Intensity", Trans. Arner. Geophys . Union, Vol. 22, 
pp. 709-721. 

Lipes , R.G. , [1979], "Altimeter Path Correction", JPL IOM 331- 
78-166A - Jet Propulsion Laboratory, Calif. Institute of 
Technology, Pasadena, California. 


Page 65 


References (Contd.) 

Lorell, J., [1979], "Altimeter Signal Anomalies Due to Sea-Surface, 
Waveheight and Tracking System Effects", Seasat Colloquim, 
October 29 - November 1, Scripps Institution of Oceanography, 

La Jolla, California. 

Lopez, R.E. , [1977], "The Lognormal Distribution and Cumulus Cloud 
Populations", Monthly Weather Review, Vol. 105, July, pp. 
865-872. 

Nagler, R.G. and D.N. Brown, [1979], "Ocean Current Monitoring 

from Space", System Planning Corporation, SPC Log No. 80-0148 
(1500 Wilson Boulevard, Suite 1500, Arlington, Virginia 22209). 

Nathanson, F.E. and P.L. Smith, Jr., [1972], "A Modified Coefficient 
for the Weather Radar Equation", in Proc. 15th Radar Meteorol. 
Conf . , Urbana- Champaign, IL, American Meteorological Society, 
pp. 228-230, October 10-12. 

Nathanson, F.E., [1969], "Radar Design Principles", McGraw Hill 
Book Co., pp. 206-209. 

Oguchi, T. , [1973], "Attenuation and Phase Rotation of Radio Waves 
Due to Rain Calculations at 19.3 and 34.8 GHz", Radio Science, 
Vol. 8, No. 1, January, pp. 31-38. 

Oh, In Hwan, [1974], "The Average Solar Flux (10.7 cm) Model For an 
11-Year Cycle", Wolf Research and Development Corporation, 

NAS 5-11933 MOD. 13, September. 

Olsen, R.L., D.V. Rogers and D.B. Hodge, [1978], "The aR b Relation 
in the Calculation of Rain Attenuation", IEEE Trans. Anten. 
and Prop., Vol. AP-26, No. 2, March, pp. 318-329. 

Ray, P.S., [1972], "Broadband Complex Refractive Indices of Ice 
and Water", Appl. Opt., Vol. II, pp . 1836-1844. 

Rogers, D.V. and R.L. Olsen, [1975], "Delay and Its Relation to 

Attenuation in Microwave Propagation Through Rain", USNC/URSI 
Annual Meeting, Boulder, Colorado, October 20-23. 


Page 6b 


References (Contd.) 

Saastamoinem, J., [1972], "Atmospheric Correction for the Tropo- 
sphere and Stratosphere in Radio Ranging of Satellites", 
Geophysical Mongraph # 15, American Geophysical Union, 
Washington, D.C. , William Byrd Press, Richmond, Va. 

Seasat (Gulf of Alaska) Workshop Report, [19791, Vol. 1, Panel 
Reports, National Aeronautics and Space Administration, 
(Published by Jet Propulsion Laboratory, California In- 
stitute of Technology, Pasadena, California, April) (G22-101). 

Seasat (Interim Geophysical Data Record) User's Handbook - Initial 
Version - Altimeter, [1979], National Aeronautics and Space 
Administration, Published by Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena, California 
91103, April, (622-97). 

Setzer, D.E., [1970], "Computed Transmission Through Rain at 

Microwave and Visible Frequencies", Bell System Technical 
Journal, Vol. 49, No. 8, October. 

Skolnik, M.I., [1974], "The Application of Satellite Radar for 

the Detection of Precipitation", NRL Memorandum Report 2896, 
October (N.R.L., Washington, D.C.). 

Smith, E.K. and S. Weintraub, [1953], "The Constants in the Equa- 
tion for Atmospheric Refractive Index at Radio Frequencies", 
Proc. IRE, Vol. 41, August, pp. 1035-1037. 

Soicher, H. , [1977], "Ionospheric and Plasmaspheric Effects in 
Satellite Navigation Systems", IEEE Trans, on Anten. and 
Prop., Vol. AP-25 , No. 5, September, pp. 705-708. 

Soicher, H. , [1975], "Plasmaspheric Contribution to Group-Path- 
Delay of Trans-ionospheric Satellite Navigation-Signals", 
AGARD Conference Proceedings No. 173, on Radio Systems and 
the Ionosphere, EM Wave Prop. Panel, Athens, Greece, 26-30 
May, pp. 6-1 to 6-15. 


Page 67 


References (Contd.) 

Tapley, B.D. , B.E. Schutz, J.G. March, W.F. Townsend, and G.H. 

Barn Editors [1979], "Accuracy Assessment of the Seasat 
Orbit and Height Measurement", Institute for Advanced Study 
in Orbital Mechanics (IASOM TR 79-5) , The University of 
Texas at Austin, Austin, Texas, October. 

Townsend, W.F. , [1980], "On Initial Assessment of the Performance 
Achieved by the Seasat-1 Radar Altimeter", NASA Technical 
Memorandum 73279, NASA Wallops Flight Center, Wallops Island, 
Virginia 23337. 

Van De Hulst, H.C., [1957], "Light Scattering by Small Particles", 
New York; John Wiley Sons, Inc., 1957. 

Wu, S.C., [1977], "Ionospheric Calibration for Seasat Altimeter", 
JPL Engineering Memorandum 315-34, Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena, California. 


